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ABSTRACT 

In the core accretion model, gas giant formation is a race between growth and migration; for a core 
to become a jovian planet, it must accrete its envelope before it spirals into the host star. We use a 
multizone numerical model to extend our previous investigation of the "window of opportunity" for gas 
giant formation within a disk. When the collision cross-section enhancement due to core atmospheres 
is taken into account, we find that a broad range of protoplanetary disks posses such a window. 
Subject headings: protoplanetary disks, origin: solar system — planets: formation, plancts:extrasolar 



1. INTRODUCTION 

Core accretion and direct gravitational instability are 
the two main competing models for how gas giant planets 
form. The core accretion picture appears to be favoured 
by both observations (Fischer & Valenti 2005) and the- 
ory (Rafikov 2005) but its viability in the face of puta- 
tive rapid orbital decay of protoplancts by planet-disk 
interactions (Ward 1997; Papaloizou & Larwood 2000; 
Tanaka et al. 2002) has long been a cause of worry. 
Thommes & Murray (2006, herafter TM06) showed that 
in fact, even with rapid migration operating, a window of 
opportunity for successful core accretion can open later 
in a protostcllar disk's lifetime, once it has dissipated to 
the point that accretion and migration timescales become 
comparable. Here, we build on this work, adopting a 
more sophisticated model which explicitly tracks the his- 
tory of each individual protoplanet all through the era of 
core formation. In this way, we are able to assess the role 
of protoplanet mergers (as opposed to just planctcsimal 
accretion) in the formation of cores, as well as to make 
predictions about the actual number of potential cores a 
given system will form. We go on to include the effect 
of gas atmospheres on the protoplanets, which increases 
their accretion cross-sections. In agreement with previ- 
ous results, we find this greatly speeds up accretion once 
protoplanet masses reach ~ 10 _1 Earth masses (M®). 
As a result, core accretion is successful, notwithstanding 
type I migration, for a wide range of disk parameters. In 
§2, we describe the model. In §3, we perform a sample 
computation using the model. In §4, we add core at- 
mospheres to the model and then recompute the sample 
case. In §5, we examine core formation likelihood as a 
function of disk properties. Conclusions are presented in 
§6. 

2. THE MODEL 

As in TM06, we adopt for the type I migration rate 
the result obtained by Tanaka et al. (2002). Likewise, 
we again use a "particle in a box" estimate of the pro- 
toplanet accretion rate, assuming it to be in the "oli- 
garchic" regime of growth (Kokubo & Ida 1998, 2000; 
Thommes et al. 2003). Also, we use the simplifying as- 
sumption of a uniform planetesimal mass m. The accre- 
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tion rate of a protoplanet is 
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where E m is the surface density of the planetesimal disk, 
H m is its scale height, Rm is the protoplanet radius, 
v csc is the escape velocity from the protoplanct's surface, 
and v re i is the average relative velocity between the pro- 
toplanet and a planetesimal. The planetesimal velocities 
are increased through gravitational stirring by the proto- 
planets, and reduced by aerodynamic drag from the gas 
disk; we take the corresponding rates from Ida & Makino 
(1993) and Adachi et al. (1976) respectively. Unlike in 
TM06, we explicitly calculate the velocity evolution of 
the plantcsimal disk, without assuming cqulibrium val- 
ues for the random velocities of the plantesimals. 

We compute the evolution of the plantcsimal disk us- 
ing a multizone approach, dividing it up into radial bins. 
Each timestep, we compute the accretion zone, as well 
as the stirred zone, of each protoplanet, noting the plan- 
etesimal disk annuli it encompasses. The stirred zone 
has a radial width 
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(Ida & Makino 1993) where tm is the protoplanet's or- 
bital radius, ru = (M/3M*) 1/,3 rjvf its Hill radius, and 
M* the mass of the central star. The accretion zone 



has a width Ar 
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of growth, adjacent protoplancts maintain a separation 
of ~ 10r#. We start the protoplanets with this spac- 
ing, and merge a pair of neighbors whenever they come 
within 7rjj of each other. At early times, such mergers 
represent the fact that there are winners and losers as 
neighboring protoplanets compete for building material. 
Caution is required in interpreting any mergers which 
occur at late times; we consider this issue in §3 below. 

3. RESULTS 

We begin by re-calculating the "baseline" model of 
TM06 using the above multizone approach. That case 
begins with a disk of mass M<j = 0.15M Q around a Solar- 
mass star. The disk has an a-parametcrized viscosity 
with a = 10 -2 , an initial outer radius of 50 AU, and a 
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Fig. 1. — Protoplanet migration (top) and growth (bottom) in 
an evolving protostellar disk with the parameters given in §3. The 
vertical line indicates f30AU; the time at which the gas disk mass 
in the inner 30 AU drops to 1 Mj up . Filled circles show when a 
core first reaches M rwy , the runaway gas accretion mass. 

scale height profile H = 0.047(r/lAU) 5 / 4 AU. The gas 
disk's time evolution is calculated analytically using the 
similarity solution of Lyndcn-Bell & Pringlc (1974). At 
t = 0, the accompanying solids disk follows the surface 
density profile of the gas, with a metallicity [Fe/H]=0.25. 
It consists of planetesimals with a uniform radius of 1 km. 
We take the "snow line" to be at 2.7 AU as in the model 
of Hayashi (1981) (sec TM06 for details). As in TM06, 
we stop the calculation when the gas mass interior to 100 
AU drops below one Jupiter mass (tiooAU = 5.5 Myrs), 
noting also the time at which the gas mass interior to 30 
AU drops below this value (t^oAU = 2.4 Myrs). As a new 
addition, we calculate the critical core mass Af cr ; t for run- 
away gas accretion as in Ikoma ct al. (2000), assuming an 
envelope grain opacity of k = 0.4cm 2 g _1 (see §4 below). 
This expression for M cr it becomes zero when planctcsi- 
mal accretion ceases, so we choose the mass for runaway 
gas accretion onset to be M rwy = min(M a it, 10M®); en- 
velope growth time for a 10 M® body not accreting plan- 
etesimals is only ~ 10 5 yrs (Ikoma et al. 2000). 

The result of the calculation is shown in Fig. 1. This 
can be directly compared to Fig. 6 of TM06. The overall 
evolution matches our previous semianalytic result quite 
closely: By t 30 A\j, the largest protoplanet mass is just 
under 10 M®, and by iiooAU it is a bit over. No pro- 
toplanets have reached M rwy by i30AU, while two reach 
this threshold between t30AU an d iiooAUi suggesting that 
core accretion growth of gas giants does have a chance to 
occur before the gas is lost. The multizone approach now 
allows us to see the role which mergers between neighbors 
plays in the growth of the protoplancts. Both the r vs. t 
and M vs. t panels of Fig. 1 clearly show these mergers. 
Because we start the protoplancts with a uniform sepa- 
ration in rn, and adopt a uniform separation at which 
to merge neighbors, the mergers occur in successive, or- 
derly generations which propagate outward through the 
disk over time. From our chosen initial protoplanet mass 
of 10~ 3 M®, three generations of mergers occur, with 
the last one taking protoplanet masses from ps 0.2 to 0.4 
M®. Thus, protoplanet mergers play a significant role in 
growth (for a discussion, see Chambers 2006), but only 
up to a mass < 1 M® in the case shown here. After that, 
the rate at which neighboring protoplanets diverge by 
type I migration exceeds the rate at which mass accre- 



tion expands their Hill radii, so that subsequent growth 
proceeds only by the sweep-up of the remaining planetes- 
imals. 

An exception to the orderly sequences of mergers takes 
place within the first 2 x 10 5 years, beginning at the snow 
line. Here, the wave of mergers starts to go inward with 
time rather than outward, with the result that one pro- 
toplanet accretes all its inner neighbors before falling off 
the inner edge of the disk. It is able to do this because the 
density jump at the snow line is steep enough that within 
it, the growth timescale decreases outward. This sets the 
scene for a "predator" protoplanet to arise which grows 
faster than its inner neighbors. Since the type I migra- 
tion rate of a body is proportional to its mass, this body 
then overtakes its neighbors, accreting them all one by 
one, and finally exiting the inner edge of the disk just as it 
attains Af rwy . However, because the growth of the preda- 
tor becomes dominated by mergers, and also because this 
is happening in a part of the disk where the majority of 
the planetesimal mass has already been swept up, we are 
no longer in the true oligarchic growth regime. Thus our 
treatment of protoplanet mergers, being statistical in na- 
ture, loses its validity. All we can say is that convergent 
migration of protoplancts takes place from the snow line 
inward. Previous N-body work suggests that differential 
migration tends to result in the establishment of mean- 
motion resonances among protoplanets (Thommes 2005; 
McNeil et al. 2005; Terquem & Papaloizou 2006) 

4. ATMOSPHERE-ENHANCED ACCRETION 

In the core-accretion model, a core becomes a gas gi- 
ant when it undergoes runaway accretion from the sur- 
rounding gas disk, thus obtaining a massive atmosphere 
(Cameron 1978; Mizuno et al. 1978). Well before that 
time, however, the core will already start to accrete some 
gas (Pollack et al. 1996). Planetesimals which intersect 
the atmosphere lose energy to aerodynamic drag; if they 
lose enough energy they are captured, even if their orig- 
inal trajectory would have caused them to miss the core 
(Pollack et al. 1996; Inaba & Ikoma 2003; Inaba et al. 
2003). Atmospheres thus provide an enhancement to 
the cores' planetesimal capture cross-sections, which in- 
creases with the density of the atmosphere and decreases 
with plantesimal size. As shown by the above authors, 
atmospheric capture can substantially boost a proto- 
planet 's accretion rate beyond the early stages of its 
growth, once it is > 10 _1 M®. 

In order to add this effect to our calculation, we use 
the capture model of Chambers (2006) which is in turn 
a simplified version of the the model of Inaba & Ikoma 
(2003). The atmosphere-enhanced capture radius Rc is 
then given by 

(Rc_ \ 3 _ 0.0790/i 4 cE M r g ( M \ 2 ( 24 \ (dM\~ 1 
\Rm) ~ nr m \M Q ) 1,24 + 51/ \dt ) 

(3) 

where \x is the mean molecular weight of the gas (taken as 
2.8, Solar composition), c is the speed of light, r m is the 
planetesimal radius, n is the atmospheric opacity, and 
e = erm/rn- We substitute Rc for Rm in our calcula- 
tion, subject to the condition that Rm < Rc < ^atmos, 
where i? a tmos is the outer radius of the protoplanet's at- 
mosphere, taken to be the smaller of its Hill and Bondi 
radii {R B = GM/c 2 s ). 
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Fig. 2. — Protoplanct migration (top) and accretion (bottom) for 
the same parameters as in Fig. 1, but including the enhancement 
to the planctcsimal capture cross-sections due to the growth of 
atmospheres on the protoplanets. 

We repeat the calculation of Fig. 1, this time with 
atmosphere-enhanced capture radii for the protoplan- 
ets. We use the opacity from the analytic fit of 
Inaba & Ikoma (2003), with a grain depletion factor of 
/ = 0.1 relative to the interstellar value to account for 
grain growth and settling (Podolak 2003). For simplicity, 
we use a single value for the opacity. We conservatively 
choose the highest (i.e. lowest-temperature) value out 
of the range: k = 4/ cm 2 g _1 =0.4 cm 2 g _1 . Fig. 2 
shows the result. Until protoplanet masses reach ~ 10 _1 
Mg, accretion proceeds in the same manner as before. 
However, once the atmospheres become effective, growth 
quickly becomes much faster than in the unenhanced 
case. Discounting the initial burst of growth at the edge 
of the snow line (see §3 above), the first protoplanet to 



M rwy does so at 4.5 x 10 5 yrs already. By t 



30AU 



2.4 



Myrs, a total of eight protoplanets have managed to 
reach Af rwy , each one as its planctcsimal accretion ceases 
and M cr ; t decreases. Since the later ones originate from 
further out in the disk, they survive for longer before mi- 
gration removes them; the last four protoplanets reaching 
M rwy before t30AU last another ~ 0.5 — 1 Myrs. Thus the 
window of opportunity for a core to become a gas giant is 
considerably larger than in the calculation without atmo- 
spheres, §3 above. It should be noted that the calculation 
does not take into account various effects which become 
important at large core mass, in particular gap-opening 
and the extra mass contributed by the atmosphere, so 
any growth which happens at M ^> 10 M ffl should not 
be taken too literally. Likewise, the series of mergers 
starting inside 15 AU at t w 4 Myrs is subject to the 
caveat discussed in §3. 

Calculations of concurrent core formation and migra- 
tion were also performed by Alibert et al. (2005), who 
concluded that type I migration would have to be re- 
duced by an order of magnitude relative to the result of 
Tanaka ct al. (2002) in order for core accretion to suc- 
ceed. The disparity between their results and ours ap- 
pears to stem from two main sources: First, they as- 
sumed a characteristic planetesimal size of 100 km rather 
than 1 km; this slows accretion relative to our calculation 
by reducing both the damping of planetesimal random 
velocities and atmosphere enhancement of protoplanet 
accretion radii. It is worth noting that in a direct com- 
parison using the same parameters, we actually obtain 




Fig. 3. — The length of the time window f w ; n (contours, in Myrs) 
for the formation of gas giant planets as a function of initial disk 
mass a viscosity parameter, and metallicity [Fe/H], The larger 
t w i n , the more likely it is that one or more gas giants will form in 
a given system; for t wln < 10 yrs, there is unlikely to be sufficient 
time. 

slower growth than Alibert et al. (2005), likely because 
their code incorporates weaker protoplanet-planetesimal 
gravitational stirring. Secondly, they chose a rather large 
initial mass of 0.6 M ffi for their protoplanets, with the 
result that at t = 0, all of them arc already migrating 
quite rapidly; throughout much of the disk, this removes 
the opportunity for disk dispersal to equalize growth and 
migration timescales. 

5. GIANT PLANET FORMATION AS A FUNCTION OF 
DISK PROPERTIES 

As in TM06, we now build up a global picture of how 
the likelihood of gas giant formation by core accretion 
changes as a function of disk properties, this time includ- 
ing the effect of core atmospheres. We compute a grid of 
models, 0.01 M < M d < 0.4 M and 10~ 4 < a < 10 -1 , 
and repeat this for metallicities of [Fc/H]=-0.25, and 
0.25. In TM06, the largest protoplanct mass at tiooAU 
was plotted for each model. Here, we instead plot the 
length of the time window, i win , between the time the 
first protoplanet reaches M rwy , and i30Au; we consider 
this to be a better measure of how likely a given system is 
to succeed in forming a gas giant. Also, we cut all compu- 
tations off at 10 Myrs, even if £3oau > 10 Myrs, since this 
is the observational upper age limit for gas disks. Fig. 
3 shows the results, indicating that t w ; n increases with 
[Fe/H], while the "optimal" value of a, which allows suc- 
cessful core accretion at the lowest disk mass, decreases 
with [Fe/H]. Within the range of obscrvationally-infcrrcd 
disk viscosities, 10~ 3 < a < 10~ 2 (Hartmann et al. 
1998), [Fe/H]=0.25 disks as low as ~ 0.03 M© in ini- 
tial mass have i w i n > 10 5 yrs, and thus stand a chance 
of accreting a Jovian envelope in time (Hubickyj et al. 
2005). For disks with [Fc/H]=-0.25, M d > 0.08 M Q is 
required. 

6. CONCLUSIONS 

We have presented here a multizone calculation of con- 
current planet growth and accretion in an evolving pro- 
toplanetary disk, focusing on whether giant planet cores 
can form in time. The first part of our results agrees well 
with our previous semi-analytic model, while elucidating 
the role played by protoplanet-protoplanet mergers at 
early times. Since the calculation tracks individual pro- 
toplanets, we can now estimate the total number of bod- 
ies which reach core mass in a given system. An impor- 
tant point is that although migration impedes planetes- 
imal accretion by moving protoplanets to the already- 
depleted inner disk before they can finish their growth, 
this actually helps in attaining runaway gas accretion, 
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provided a protoplanet is able to reach > 10 M s before 
its planetesimal supply is cut off. 

When we add to our calculation a simple model for 
the effect of gas atmospheres on growing protoplanets, 
the results change dramatically. In our baseline model, 
the time to grow a body large enough to serve as a 
giant planet core shrinks from over 3 million years to 
less than half a million years. Thus in the competi- 
tion between accretion and migration, the balance swings 
strongly toward the former, relatively speaking. This 
makes the window of opportunity for the formation of gas 
giants by core accretion — the time between the appear- 
ance of the first core-sized body, and the disappearance 
of the gas disk — much larger for a given set of parame- 
ters. As a result, the range of disk properties which can 
plausibly produce gas giants, notwithstanding effective 
type I migration, becomes significantly broader than in 



the atmosphere-less calculation of TM06. Now, we find 
that a Solar-metallicity disk of moderate initial mass 
(Md > 0.05 M Q ) stands a reasonable chance of giving 
birth to a jovian planet. Since we neglect fragmentation 
of planctcsimals, and the potentially very rapid accre- 
tion of fragments at low relative velocities dominated by 
Kcplcrian shear (Rafikov 2004), which will also further 
increase the role of atmospheric capture (fnaba & Ikoma 
2003), the accretion rates in our calculation are likely 
conservative, and core formation may thus proceed even 
more readily; these effects will be added in future work. 
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